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2.  Emission  of  Hot  Electrons  from  Silicon  into  Silicon  Dioxide 

Our  model  calculations  of  this  effect  have  shown  that  the  emission 
takes  place  because  of  collision  broadening  effects  due  to  the  electron- 
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continuity).  Our  findings  have  been  confirmed  experimentally  by  Muller 
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MOS  transistors. 

3.  Hot  Electron  Shockley  Equations 
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Equations  and  include  in  a  straightforward  fashion  hot  electron  effects 
(also  diffusion). 
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MSS*-  -  Hoa»  **» aa  dCveloppd  uaa  aiaailacioa  par  la  mfthode  da  Mont*  Carlo  pour  du 
silieium  at  da  1  AaCa  as  ineluaae  uaa  * true Cura  da  banda  rdaliata.  Laa  taux  d'ioni- 
aacioa  par  impact  at  laa  vital aaa  da  ddriva  aa  rdgime  coacinu  aoua  forts  chaapa  (lac- 
criques  (>100  kV/ca)  oat  ltd  calculds  k  diff Creates  taapdraturaa. 


Aba tract.  -  Wa  hava  performed  a  Mont a  Carlo  alaulatlon  for  GaAa  and  SI  with  the 
raaliatlc  band  atructura  includad.  Staady  atata  Impact  Ionization  rates  and  drift 
velocities  under  high  electric  fields  (>  100  kV/ca)  ware  calculated  at  various 
temperatures. 

1.  Introduction.  -  Impact  ionization  is  an  essential  mechanism  In  the  operation  of 
semiconductor  devices  such  as  avalanche  photo-diodes  or  transit  flaw  devices.  The 
dependence  of  Impact  Ionization  on  the  crystallographic  orientation  has  attracted 
substantial  Interest  because  of  its  relevance  to  noise  and  other  phenomena  in  these 
devices.  This  dependence,  however,  ia  not  shown  by  any  of  the  theories  aa  given  by 
Wolff  (1],  Shockley  (21,  andiBaraff  [3],  since  none  of  then  include  a  realistic 
bead  structure. 

Ws  have  developed  a  complete  theory  for  lapact  ionization  and  generally  high 
field  transport  la  semiconductors  by  containing  a  Nonce  Carlo  simulation  with  the 
realistic  bead  structure  calculated  by  the  enpirleal  pseudopotential  method.  Ws  do 
take  Into  account  scattering  by  all  possible  phonon  types,  the  change  in  the  density 
of  states  high  In  the  band,  the  exact  velocity  v  ■  ^  V^E(S)  (no  effective  mass 
approximation),  the  collision  broadening  of  the  electronic  states,  and  the  tempera¬ 
ture  effect. 

Details  of  the  nodal  and  the  results  for  GaAa  at  300  K  con  be  found  In  two  of 
our  previous  papers  [4,3].  In  tha  case  of  SI,  the  first  two  conduction  bands  were 
included,  lesldes  the  X-X  scattering  we  also  include  the  X-L  scattering  In  SI. 

In  Section  2,  we  describe  briefly  our  model  and  point  out  Its  differences  from  the 
comnonly  used  model.  The  effects  of  the  Inclusions  of  the  second  band  end  the 
transition  from  X-L  In  SI  are  discussed  In  Section  4. 

2.  Theoretical  Model.  -  The  model  for  Monte  Carlo  simulation  has  two  main  Ingredi¬ 
ents:  (1)  the  band  structure  and  (11)  the  scattering  rate.  We  describe  briefly  in 
the  following  the  different  features  that  have  been  Included  In  our  model  and  the 
advantages  It  has  over  other  models. 
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We  present  a  theory  of  deformation  potential  electron-phonon  scattering  in  the  bulk,  in 
superiattices,  and  at  interfaces  of  semiconductors.  Expressions  for  the  acoustical  deformation 
potential  constants  are  obtained  in  closed  form  for  bulk  semiconductors  and  for  a  monolayer 
superiattke.  Deformation  potential  constants  are  evaluated  numerically  for  superlattices  of 
various  thicknesses.  The  results  are  in  good  agreement  with  available  experiments  and  provide  a 
detailed  understanding  of  deformation  potential  scattering. 

PACS  numbers:  72.10.  —  d,  73.40.  —  c,  72.10.Di 


One  of  the  dominant  scattering  mechanisms  of  carriers  in 
semiconductors  involves  long  wavelength  acoustical  phonons. 
It  was  first  pointed  out  by  Bardeen  and  Shockley1  that  such 
interactions  can  be  described  by  a  deformation  potential.  The 
basic  assumption  b  that  the  local  lattice  deformation  produced 
by  the  phonons  is  equivalent  to  a  homogeneous  deformation 
of  the  crystals.  The  strength  of  this  coupling,  the  deformation 
potential  constant  (Z„(k)]  is  then  determined  by  the  shifts  in 
the  (nth)  energy  bend  (at  momentum  k)  per  unit  of  dilation 
({  «  £V/V)  produced  by  the  acoustical  waves, 

u> 

An  understanding  of  deformation  potential  scattering  is.of 
fundamental  importance  for  transport,  optical,  and  other 
properties  of  bulk*  semiconductors  and  superiattices.  Despite 
advances?  in  electronic  band  structure  calcdrtions,  there  lave 
been  only  a  few  calculations4  of  the  deformation  potential 
constants  for  a  few  selected  bulk  semiconductors  and  nooe  in 
superiattices.  The  main  difficulty  is  that  the  required  absolute 
shtft  of  the  band  edge  is  beyond  most  empirical  band  structure 
calculations  where  the  relative  band  edges  are  fitted  to  agree 
with  experimental  observations.  (See  Kleinman,  Ref.  4.)  As 
a  consequence,  the  deformation  potential  constants  for  the 
bulk  are  usually  deduced  from  transport  and  optical  mea¬ 
surements  with  various  fitting  schemes,  often  resulting  in  a 
range  of  experimental  estimates.3 

In  this  paper  we  sumiAarize  a  theory  for  the  deformation 
potential  scattering  by  long  wavelength  acoustical  phonons 
in  bulk  semiconductors3;  and  we  extend  the  theory  to  super- 
lattices  and  interfaces.  We  show  that  the  absolute  shift  of  the 
bund  edge  can  be  obtained  using  an  equilibrium  sum  rule. 
Expressions  for  the  deformation  potential  constants  are  ob¬ 
tained  in  doeed  form  for  bulk  semiconductors  and  the  results 


are  in  very  good  agreement  with  experimental  results.  The 
theory  is  applied  to  semiconductor  superiattices  and  closed 
form  expressions  for  the  deformation  potential  constants  are 
given  for  monolayer  superiattices.  Results  for  lattice-matched 
AlAs/GaAs  superiattices  of  various  thickness  are  obtained 
numerically.  To  our  knowledge  this  is  the  first  calculation  of 
deformation  potential  constants  for  interfaces  and  superiat¬ 
tices  using  realistic  models  of  band  structure.  We  show  that 
for  most  practical  cases  the  differences  between  superlattice 
and  bulk  deformation  potential  constants  are  small.  The  dif¬ 
ferences  that  do  occur  include  a  new  "edge”  scattering 
mechanism  and  are  caused  by  bond-edge  discontinuities. 

The  electronic  energy  bands  of  a  system  of  ions  and  elec¬ 
trons  can  be  written  as, 

E.(k)-E.*»(k)  +  E*  (2) 

where  E°  is  related  to  the  ion-ion  self -energy  and  is  usually 
chosen  such  that  the  zero  of  E*  “(k )  is  at  the  maximum  of  the 
valence  bond.  E0  contains  contributions  involving  ions  and 
electrons  (see  Kleinman,  Ref.  5).  It  is  density  dependent  and 
is  very  difficult  to  evaluate  in  general.  Under  dilation  or 
changes  in  the  lattice  constants,  a  shift  of  the  energy  bands  is 
produced, 

6En(k)  •  6En's(V)  +  6E°.  (3) 

Although  £°  is  difficult  to  calculate,  small  changes  in  E° 
can  be  easily  approximated  by  an  extension  of  techniques  used 
by  Chadi7  in  his  ground  state  calculations  involving  total- 
energy  minimization  of  surface  atomic  structures.  One  con¬ 
siders  a  sum  over  the  occupied  single  particle  states. 


Etot  ■  E  EBK(k)  +  NME0.  (4) 

ft.k 
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A  theory  of  deformation  potential  electron— phonon  scattering  coeffi¬ 
cients  is  presented  for  elemental  and  compound  semiconductors. 
Expressions  for  the  acoustical  deformation  potential  constants  at 
symmetry  points  are  obtained  in  closed  form  for  direct -gap  bulk  semi¬ 
conductors.  The  deformation  potential  vs  k  is  predicted  for  GaAs. 


IN  THIS  PAPER  we  present  a  simple  but  successful 
theory  of  dilations!  deformation  potentials  for  semi¬ 
conductors,  obtaining  analytic  expressions  for  bulk 
deformation  potential  constants  in  terms  of  empirical 
tight-binding  energy  band  parameters. 

The  idea  of  a  deformation  potential  was  first 
enunciated  by  Bardeen  and  Schockky  J1  ],  who  pro¬ 
posed  that  kmg-waveiength  acoustic  phonons  be  simu¬ 
lated  by  homogeneous  lattice  deformations,  permitting 
one  to  calculate  the  electron— phonon  scattering 
minx  cmamu  or  cwscviMQOa  potinttii  consuMt 
Z,(k)  by  expanding  the  band  structure  to  first  order  in 
Urn  change  et  lattice  constant  dc£: 


£<k,efc  +  dft)  -  £;(k,«i)+Z*(kX3del/«x.). 


Thus  ths  deformation  potential  conrtant  Z,(k)  is  pro¬ 
portional  to  tht  derivative  of  the  energy  band  with 
respect  to  the  lattice  constant  e*,.  (Here,  for  clarity  of 
presentation,  we  have  avoided  the  tensor  character  of 
2^(k)  by  cottskkriag  only  purely  dilations!  strains: 

|  tedf7 V*  34*1/01;  e  complete  dtsaierion  of  the 
tensor  character  wiB  be  published  separately.) 

As  riwwn  by  Kktesnan  (2],  tha  one  main  difficulty 
hi  computing  deformation  potential  constants  arises 
b  scenes  of  a  difltcnit-fo-evahMte  many-body  seif-energy 
contribution  (3],  We  circumvent  the  direct  calculation 
of  the  many-body  term  by  usings  sum  rule  (3, 4] ' 
expceesit  in  terns  of  easily  calculable  one-body 


To  iBustrate  this  point,  wt  consider  the  simple 


example  of  the  following  one  electron  Hamiltonian: 

H  -  p]pm  +  V(ri;  R)  +  S(r,;  R). 

Hen  the  crystal  potential  is  V(r;  R),  the  ion  coordinates 
are  denoted  R,  and  the  self-interactions  are  contained  in 
the  self-energy  Sfa;  R).  The  band  structure  fR(k)  can  be 
calculated  for  this  model, 

{p*/2m+  V(r;R)  +  S(tl;R)}f+(t)  -  £«(k)^(r). 

or,  as  in  the  present  case,  determined  by  empirically 
fitting  data.  Since  SQStt ;  R)  is  unknown,  £*( k)  is  not 
fully  specified.  Only  e„(k)  ■*  £„(k>— 5*  is  determined 
by  the  usual  bare  band  structure  calculation.  Homo¬ 
geneous  deformation  of  foe  lattice  changes  both  the 
crystal  potential  V  snd  the  self-energy  S.  However,  the 
changes  of  foe  sverage  self-energy  resulting  from  defor¬ 
mation  cm  be  determined  by  the  requirement  that  the 
total  energy  be  a  minimum  at  equilibrium,  as  demon¬ 
strated  by  Chad!  [4]. 
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whe-  ie  primed  sums  are  restricted  to  the  N  occupied 
H  1  states  snd  we  have  made  a  crude  approximation 
to  foe  total  energy.  Thus  the  deformation  potential 
constant  is 
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Predictions  of  the  densities  of  states  an  given  for  various  forms  of  controlled  random 
disorder  in  artificial  multilayered  materials  or  superstructures. 


L  INTRODUCTION 

With  the  development  of  techniques'  for  fabri¬ 
cating  artificial,  multilayered,  periodic  materials,  it 
is  passible  to  study  controlled  randomness  and  its 
effects  on  electronic  states  of  matter.  A  hint  of 
the  suitability  of  such  superstructured  materials  for 
randomness  studies  was  contained  in  the  early 
work  of  Esaki  and  Tsu,2  but  since  then  the  greatest 
experimental  effort  has  been  to  remove  randomness 
from  superstructures  and  to  make  dm  interfaces 
between  layers  abrupt  and  atomically  smooth. 

Superstructures  with  two  types  of  randomness 
can  be  grown:  randomly  varying  layer  thicknesses 
and  random  layer  compositions.  Both  metallic  and 
semiconductor  superstructures  have  been  grown, 
and  the  ideas  of  this  paper  will  be  applicable  to  ei¬ 
ther  however,  for  definiteness  we  confine  our  at¬ 
tention  to  superstructures  composed  of  layered  - 
semiconductors  GaAs-AlzGa,  As,  where  the  ran¬ 
dom  variables  are  t be  layer  thicknesses  d  and  alloy 
compositions  x.  We  envision  that  the  stochastic 
variables  d  and  x  are  determined  by  a  random 
number  generator  during  the  superstructure  growth 
process  and  that  these  quantities  are  preserved  for 
the  analyses  of  data  taken  from  the  superstructured 
sample. 

The  resulting  disorder  is  controlled,  finite  in  ex¬ 
tent,  and  essentially  one  dimensional.  Its  con¬ 
trolled  nature  is  a  valuable  aid  to  understanding 
random  systems,  because  few  such  systems  have 
been  fabricated  before.  Moreover,  controlled- 
disorder  superiattices  offer  the  possibility  of  study¬ 
ing  heretofore  unimagined  combinations  of  order 
and  disorder,  such  as  periodically  stacked  alternat¬ 
ing  ordered  and  disordered  arrays.  The  finiteness 
of  the  artificially  produced  disorder  offers  oppor- 
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'  (unities  to  directly  determine  the  extent  of  iocall 
tion  of  states  and  to  study  the  impending  onset  1 
Anderson  localization3  as  a  function  of  increasim 
size.  One  dimensional  random  systems  are  com| 
paratively  well  understood  theoretically;  numerolj 
models  have  been  solved  exactly  and  even  more  I 
have  been  thoroughly  studied.  I 

For  electronic  states  in  semiconductor  super-  1 
structures  three  theoretical  regimes  present  thetn-1 
selves:  the  regime  of  localized  deep-trap-lilce  1 
states,4  the  mixed  regime  of  localized  and  extends^ 
states,  and  the  regime  of  extended  states  described 
by  effective-mass  theory.3  The  most  interesting- isl 
the  one-band  effective-mass  theory  of  carrier  1 
motion,  which  we  shall  consider.  This  regime  has 
two  subregimes:  the  quantum  well  limit6  in  which 
the  effective-mas*  electron’s  de  Broglie  wavelength 
is  comparable  with  or  larger  than  typical  super-  . 
structure  dimensions,  and  the  classical  limit  in 
which  the  electronic  spectrum  is  characteristic  of  a 
classical  particle  colliding  with  barriers.  In  this  : 
paper  we  shall  restrict  our  attention  to  the  qusn-  , 
turn  well  limit”  - 

The  standard  theoretical  questions  to  be  ) 

answered  by  any  theory  of  disorder  are  as  follows: 
(i)  What  is  the  ensemble-averaged  density  of  states 
for  the  random  system?  <ii)  Are  the  band  gaps  of 
the  ordered  structure  preserved  in  spite  of  the  dis¬ 
order  or  annihilated  by  it?  (iii)  Can  selected  types 
of  disorder  introduce  gaps  into  a  spectrum  that 
would  otherwise  be  continuous?  (iv)  To  what  ex¬ 
tent  does  the  disorder  produce  localized  states,  is 
diffusion  possible,  and  what  is  the  transmission 
coefficient  for  an  electron  in  a  random  superstruc¬ 
ture?  (v)  What  are  the  effects  of  “many  body”  in¬ 
teractions  between  electrons  confined  to  adjacent 
layers  and  can  they  produce  one-dimensional 
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Ahetract— We  present  a  simple  formalism  to  include  hot  electron  diffusion  into  device  modeling  programs.  The 
effectiveness  of  this  method  is  illustrated  for  a  DMOS-structure  and  the  permeable  base  transistor.  However,  the 
formalism  is  general  and  can  be  applied  to  any  semiconductor  device. 


'  l  omoovcnoN 

It  is  well  known  that  the  drift  velocity  of  electrons 
saturates  at  high  electric  fields.  This  is  a  consequence  of 
the  high  carrier  energies  and  increased  scattering  rates. 
Stratton[l]  derived  a  simple  formalism  to  include  also 
the  changes  of  the  diffusivity  O  in  high  electric  fields.  He 
showed  that  if  a  Maxwellian  distribution  with  a  small 
drift  term  is  assumed,  the  diffusivity  is 

(1) 

which  is  formally  identical  to  the  Einstein  relation.  Here 
H  is  the  mobility,  k  is  Boltzmann’s  constant  and  T,  is  the 
carrier  temperature,  not  the  temperature  of  the  crystal 
lattice.  Stratton  derived  (1)  with  the  assumption  that  the 
scattering  rate  depends  on  the  spatial  coordinate  only 
through  the  carrier  temperature.  His  assumptions  also 
exclude  anisotropies  in  die  diffusion  constant  as  well  as 
nonlocal  effects  introduced  for  example  by  time  depen¬ 
dent  friction.  These  complications  will  be  discussed  at 
the  end  of  the  paper  where  some  important  references  to 
more  exact  treatments  are  given. 

The  main  purpose  of  the  paper  is  to  present  a  method 
which  maximizes  the  physical  input  in  the  treatment  of 
electronic  diffusion  in  device  modelling,  while  simul¬ 
taneously  minimizing  the  necessary  additional  com¬ 
putational  effort 

Although  the  method  we  present  is  more  general,  it  is 
meant  for  application  to  silicon  devices  with 
dimensions  *  0.23  x  |0~4  cm  and  for  time  scales* 
3  x  I0'w  sec. 

The  key  point  to  note  in  (I)  is  that  T,  is  an  increasing 
function  of  F-J  while  n  is  a  decreasing  function  of  F •  j. 
F  is  the  electric  field  and  |  is  the  current  density.  In  many 
device  models  the  assumption  has  been  made  that  the 
carrier  temperature.  T„  is  equal  to  the  lattice  tem¬ 
perature,  Tl.  This  has  often  given  the  impression  that 
diffusion  currents  are  unimportant[2-6).  However,  for 
high  electric  fields  and  large  charge  density  gradients  this 
assumption  causes  considerable  error  in  the  calculated 


current  densities.  Below  we  demonstrate  the  effect  of 
properly  including  the  carrier  temperature  and  mobility  m 
as  function  of  j-F  and  show  the  significance  of  this 
effect  for  both  the  diffusion  current  and  the  hot  electron 
mobility  in  devices.  The  Permeable  Base  Transistor 
(PBT)  and  the  Double-Diffused  MOS  transistor  (DMOS) 
will  serve  as  examples  because  they  exhibit  both  high 
electric  fields  and  substantial  charge  gradients.  No  • 
attempt  was  made  to  precisely  determine  the  diffusion 
constant  Thus  these  analyses  are  qualitative  rather  than 
quantitative  in  nature  and  the  results  should  not  be 
considered  exact  The  advantage  of  our  model  is  that 
computation  time  in  device  modelling  programs  should 
not  greatly  increase  when  our  procedure  is  used. 

X  ANALYTICAL  CONStDOUTIONS 

Part  of  the  following  treatment  of  high  field  diffusion 
in  devices  has  been  suggested  before  for  CCD’s[7]-  We 
outline  the  approach  here  in  more  detail  and  emphasize 
the  generality.  It  is  important  to  notice  that  the  mobility 
and  diffusion  constant  are  not  single  valued  functions  of 
the  electric  field  F.  The  reason  is  that  the  carrier  energy, 
which  ultimately  determines  p  and  D.  is  a  function  of 
J  •  F,  the  power  input  per  unit  volume  as  mentioned  in  the 
introduction. 

For  carrier  concentrations  n  >  10ITcm~’  an  electron 
temperature  Tc  can  be  defined  and  the  mobility,  as  well 
as  the  diffusion  constant,  becomes  a  function  of  Tc-  For 
R<IOITcm'3  an  electron  temperature  Tr  cannot  be 
defined  as  the  carrier  distribution  function  becomes  non- 
Maxwellian.  However,  one  still  can  define  an  average 
electron  energy  and  the  considerations  below  still  apply 
except  for  the  appearance  of  statistic  factors  which  leave 
the  treatment  qualitatively  but  not  quantitatively  valid. 
Let’s  assume  for  the  moment  that  T,  can  be  defined.  The 
characteristic  functional  dependence  of  the  mobility  on 
Tr  in  silicon  is  then  (8] 

n*  IhtS/iTJT'.).  (2) 

Here,  TL  is  the  temperature  of  the  crystal  lattice;  eqn  (2) 
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Investigation  of  Transient  Electronic  Transport 
in  GaAs  Following  High  Energy  Injection 
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Abstract -W»  present  Mont*  Cado  Simula  boos  of  the  transient  be-  by  Ruch)  to  simulate  injection  over  a  (hetero)barrier.  The 


havior  of  electrons  injected  into  CaAs  at  Irish  enetfiea  and  accelerated 
(decelerated)  by  constant  electric  fields.  Our  calculations  differ  from 
previous  calculations  due  to  the  inclusion  of  a  realistic  band  structure 
(empirical  paeudopotential)  and  the  injection  of  electrons  at  Irish 
enetpiea  (e^,  via  a  heterohanier).  The  results  show  that  a  narrow 
"colHsion-free  window  "  (CFW)  exists  with  respect  to  parameters  such 
as  the  external  electric  field,  the  injection  enersy,  the  external  vot tapes, 
and  the  semiconductor  dinsenskme.  Within  this  window  avefape  elec¬ 
tron  velocities  of  ~8  x  107  cm/s  can  be  achieved  over  distances  of 
10~*  cm  in  emitter  (source)-  and  hase-like  structures.  Voltapt  (field) 
parameters  typical  for  the  collector  (drain)  am  far  outside  the  CFW  and 
allow  only  for  much  reduced  (by  a  factor  of  ~10)  electron  vaiodtiea. 
We  also  discuss  thermal  noise  in  “ballistic  devices’*  and  show  that  the 
noise  equivalent  temperature  can  be  exceedinply  Irish. 


I.  Introduction 

SINCE  THE  FIRST  calculation  of  transient  electron  trans¬ 
port  in  short-channel  FET’s  by  Ruch  [1],  overshoot  phe¬ 
nomena  of  the  electron  velocity  on  very  short  time  and  length 
scales  have  attracted  considerable  attention  [2) -[5] .  Shur 
(3]  and  Shur  and  Eastman  [S]  added  features  of  space  charge 
limited  conditions  to  the  velocity  overshoot.  They  investigated 
the  initial  transient  and  called  it  the  “near  ballistic"  regime. 
The  term  “ballistic”  is  difficult  to  define  and  is  currently 
applied  to  a  wide  range  of  device  parameters  and  dimensions. 
Investigations  in  [3]  and  [S]  indicate  that  “near  ballistic” 
transport  over  larger  distances  (>1000  A),  is  achievable,  necessi¬ 
tates  injection  of  electrons  at  higher  energy.  High  energy 
injection  was  also  discussed  by  Hess  [6]  for  low  temperature 
transport  free  of  scattering  events  over  extremely  large  distances 
L  >  10*4  cm  (for  electron  energies  below  0.036  eV). 

The  criteria  for  designing  devices  and  choosing  materials  such 
that  high  transient  speeds  can  be  advantageously  achieved  using 
high  energy  injection  are  the  key  issue  of  this  paper.  We 
approach  the  problem  by  a  Monte  Carlo  simulation  which 
includes  the  details  of  the  band  structure  as  calculated  by  the 
empirical  pseudopotential  method.  The  electron  is  started  at 
higher  energy  values  (not  at  the  bottom  of  the  band  as  done 
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electric  field  is  chosen  to  be  constant  over  the  whole  distance 
of  the  simulation,  in  contrast  to  the  choice  of  the  self-consis¬ 
tent  field  by  Shur  [3].  The  reason  for  this  is  the  following: 
Within  the  CFW,  the  electron  velocity  does  not  vary  strongly 
over  the  device  length  because  the  electrons  are  already  injected 
at  high  velocities  (in  contrast  to  the  cases  considered  by  Shur). 
Therefore,  the  carrier  concentration  is  rather  constant  and  the 
electric  field  induced  by  the  carriers  is  of  minor  importance. 
In  all  other  respects  our  calculation  is  standard.  We  do  not 
include  the  intracollisional  field  effect  [7]  because  it  would 
greatly  complicate  the  computations.  In  the  CFW  it  also 
would  make  only  minor  contributions  since  the  electric  fields 
in  this  window  are  small. 

With  respect  to  device  applications,  the  result  of  the  calcula¬ 
tions  can  be  summarized  as  follows:  On  a  length  scale  of  1000 
A,  emitter  (source)-  and  base-tike  structures  may  show  effects 
typical  for  collision-free  transport  (if  carefully  designed); 
collector(drain>like  structures  will  not. 

II.  Physical  Model  and  Method  of  Computation 

As  discussed  above,  we  consider  high  energy  electrons  injected 
into  GaAs  (e.g.,  from  AlxGa,.JtAs  or  6-like  electric  fields 
created  by  spice-charge  layers).  The  transition  is  assumed  to 
be  abrupt,  i.e.,  the  electron  gains  the  kinetic  energy  A£  and  a 
forward  momentum  AJc  when  transferring  to  the  GaAs  with¬ 
out  any  energy  losses. 

As  illustrated  below,  the  calculations  of  the  self-consistent 
field  resulting  from  external  voltages  and  electron  redistribu¬ 
tion  in  the  GaAs  is  a  difficult  problem  and  depends  on  many 
details,  most  importantly: 

i)  the  boundary  conditions  of  the  injection, 

ii)  the  statistics  of  electrons  and  impurities  in  the  GaAs, 

iii)  the  level  of  injection,  and 

iv)  the  velocity  distribution  and  the  injection  energy. 

Let  us  assume  that  our  device  is  short  in  the  x -direction  but 
rather  wide  in  the  ^-direction.  We  then  can  stiii  define  average 
quantities  such  as  the  density  of  electrons  in  a  meaningful  way 
and  use  a  continuum  picture  as  follows:  If  the  distribution 
function  is  denoted  by  /we  define  electron  density 


or  the  current  density  as 
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We  present  the  results  of  a  theory  of  electronic  energy  bands  of  super¬ 
lattices  using  an  application  of  a  recent  novel  empirical  tight  binding 
method.  The  theory  is  applied  to  GaAs-AlAs  (001)  superlattices  with  very 
good  agreements  with  available  experimental  results  for  both  the  direct 
and  indirect  energy  gaps. 


SEMICONDUCTOR  SUPERLATTICES  of  alternating 
layers  of  GaAs  and  AlAs  have  been  the  subject  of  exten¬ 
sive  investigation  [1].  These  superlattices  provide  a 
means  of  studying  properties  of  electrons  confined  to 
two  dimension  as  well  as  new  materials  with  properties 
very  much  different  than  their  built  constituents.  Some 
of  the  experimental  results  [2]  can  be  explained  on  the 
basis  of  a  one  dimensional  Kronig-Penny  model  where 
the  AlAs  layers  are  assumed  to  form  potential  barriers 
confining  the  electrons  and  holes  of  GaAs.  The  barrier 
heights  are  obtained  from  the  band  edge  discontinuity. 
For  a  more  detailed  understanding  of  the  electronic 
structure  of  these  superlattices,  two  calculations  [3, 4] 
using  empirical  pseudopotential  (EP)  and  one  calcu¬ 
lation  [5]  using  empirical  tight-binding  (ETB)  have  been 
reported.  These  calculations  have  been  recently  dis¬ 
cussed  [6]  and  compared  to  experimental  results.  The 
authors  (6]  noted  that  only  one  of  these  calculations 
[4]  gives  results  in  reasonable  agreement  with  exper¬ 
imental  data.  In  this  communication,  we  present  the 
results  of  a  tight-binding  theory  for  superlattice  using 
an  application  of  a  recent  novel  empirical  tight-binding 
method.  Results  for  the  direct  and  indirect  band  gaps 
of  a  range  of  monolayer  of  AlAs  and  m  layers  of  GaAs 
superiattices  [(GaAs)^(AlAs)|]  form  -  l,2,3,and8 
ate  compared  to  experimental  results.  The  theory  pro¬ 
vides  a  good  description  of  the  energy  gaps  over  the 
entire  range  of  these  superiattices. 

The  tight-binding  methods  as  applied  to  super- 
lattices  with  alternating  layers  of  GaAs  and  AlAs  have 
been  discussed  in  detail  by  Schulman  and  McGill  (SM) 
(5(.  Their  methods  employ  an  enlarged  unit  cell  which 
then  describes  the  superlattice  periodicity  in  the  (001 ) 
direction.  The  atoms  in  each  layer  are  described  by 
tight-binding  parameters  of  the  bulk  using  a  Spi  basis 
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with  nearest  and  second  nearest  neighbor  interactions. 
The  band  edge  discontinuity  is  incorporated  by  subtract¬ 
ing  the  valence  band  discontinuity  from  the  AlAs 
diagonal  bulk  parameters.  The  second  nearest  neighbor 
interactions  were  used  to  improve  the  shape  of  the  bulk 
conduction  band. 

We  have  used  a  method  similar  to  that  of  SM  with 
some  exceptions.  We  used  a  recent  novel  empirical  tight- 
binding  method  [7]  which  is  able  to  provide  good 
description  of  the  conduction  bands  as  well  as  the 
valence  bands  of  bulk  semiconductors  by  introducing  an 
excited  S  state  into  the  Sp*  minimum  basis  with  nearest 
neighbor  interactions.  This  method  provides  a  conveni¬ 
ent  and  reasonably  accurate  representation  of  the  bulk 
band  structure  for  a  wide  range  of  semiconductors  [7]. 
As  in  SM,  the  valence  band  edge  discontinuity  (2]  is 
taken  to  be  15%of  the  direct  gap  differences  of  GaAs 
-and  AlAs  and  is  subtracted  from  the  AlAs  diagonal  bulk 
parameters.  The  resulting  matrix  of  dimension 
(n  +  m)  x  10  [for  (GaAs^fAlAs),,  superiattices]  is 
diagonalized  numerically  to  obtain  the  electronic  energy 
band  structure.  The  available  experimental  data  on  the 
direct  and  indirect  band  gaps  are  very  few.  To  test  our 
theory,  we  have  performed  calculations  on  the  series  of 
monolayer  superiattices  (GaAs-AlAs),  for  which  very 
good  experimental  results  have  been  explicitly  given  [8] 
(for  m  ■  l,  2, 3,  and  8)  including  the  indirect  band  gap 
for  the  m  ■  1  case.  See  Fig.  1 .  Given  the  theoretical  and 
experimental  uncertainty  (such  as  accuracy  in  desenp- 
tion  of  the  bulk  band  structure,  band  edge  dir  continuity 
and  layer  thickness  fluctuation),  the  rather  gjod  agree¬ 
ment  for  the  m  ■  l  case  (both  the  measured  direct  and 
indirect  pps  are  close  to  the  calculated  valves)  may  be 
fortuitous.The  overall  agreement  over  the  wide  tangs 
of  m  ■  1 , 2, 3,  and  8,  however,  indicates  that  our  theory 
indeed  can  provide  a  good  description  of  the  super  lattice 
energy  bands  within  the  Slater -Koster  (9)  interpolating 
scheme.  As  was  noted  in  (6),  theix  EP  calculations  (4) 
differ  from  the  experimental  rerdts  by  about  60meV 
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HIGH  ENERGY  DIFFUSION  EQUATION  FOR  POLAR  SEMICONDUCTORS 


J.  P.  Leburton  and  K.  Hess 


Department  of  Electrical  Engineering  and 
•  Coordinated  Science  Laboratory 

University  of  Illinois  at  Irbana-Cliarapaign 
Urbana,  Illinois,  I’.S.A.  61S01 

:.«’e  present  a  novel  method  to  explicitly  solve  the  Boltzmann  equation  for  higiilv  ener¬ 
getic  electrons  interacting  with  polar  optical  phonons  and  scattering  mainly  in  for¬ 
ward  direction.  Our  approach  reduces  the  collision  integral  of  the  Boltzmann  equation 
Into  a  differential  operator  which  is  much  easier  to  manipulate  than  the  integral  form 
and  does  not  require  a  relaxation  time  approximation.  The  spatial  diffusion  of  highly 
energetic  electrons  is  calculated  and  discussed  In  context  with  high  speed  velocity 
transients  t"ballistie  transport").  Explicit  results  compare  favorably  with  sophisti¬ 
cated  Monte  Carlo  simulations  and  are  well  suited  to  treat  complex  transport  problems 
in  submicron  III-V  devices.  As  an  exanple  we  calculate  tite  current  voltage  character¬ 
istic  of  a  "camel-back  diode". 


In  polar  semiconductors  and  for  highly  energetic 
processes,  the  scattering  of  the  charge  carriers 
hy  polar  optical  phonons  (P.O.P.)  is  mainly  in 
the  forward  direction,  i.e.,  with  conservation 
of  the  electron  momentum  direction.  'Under  these 
circumstances,  the  collision  Integral  of  the 
Boltzmann  equation  is  reduced  ro  a  differential 
operatur  [i,2| 
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which  Is  much  easier  to  manipulate  chan  Che 
integral  form.  Here  c  is  the  carrier  energy, 
1/'  is  the  spontaneous  P.O.P.  emission  rate, 

.  is  the  P.O.P.  frequency,  and  '  is  a  factor 
depending  on  the  temperature  T,  given  by 


1 


-z 
e _ 


1  +•  e 


-* 


with  z  •  ft-/k  T. 
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Figure  I:  Spatial  evolution  distribution  at 
J«JQOK.  The  starting  distribution  function 
(•-function)  is  represented  by  the  vercical  line 
at  ■  .*  yli  ,  1,2,3,  and  -l  are  for 


0.25,2,1,  and  *  respectively. 


Solid  line:  equilibrium  distribution  given  bv 
Eq.  (3) 

Dotted  line:  Maxwell-Boltzmann  distribution 


The  differential  nature  of  the  collision  opera¬ 
tur,  which  is  of  Fokker-Planck  type,  is  well 
adapted  to  treat  the  stochastic  processes  of 
non-equilibrium  phenomena.  Moreover,  the 
solution  of  the  Boltzmann  equation  does  not 
require  the  relaxation  tine  approximation,  and 
mere  gcner :  1  solutions  can  easily  be  derived. 

A  specific  application  of  our  cheorv  is  the 
important  case  ibal llstlc  transport  )  of  steady- 
state  and  field-free  spatial  dit fusion  of 
.  aar.te  injected  at  high  energy,  in  Fig.  I,  the 
reiax.it  leu  ei  an  xenoenereec  tc  distribution 
i  UU'  tin  iD.F.)  .mproxlnated  by  a  -function  ac 
tin  “initial”  position  x  *  ■»  is  shown  for  a 
p..rabol  i,  hand.  lypicallv.  tile  evolution  of  Che 
,).i.  au  oe  described  as  tel  levs.  initial  Iv, 
tie-  is  sir ’Ugly  peaked  at  »  ,  then 


under  the  influence  of  the  P.O.P.  scattering, 
mainly  emission.  It  drifts  toward  low  energy 
and  is  broadened.  The  "drift  coefficient”  in 


energy  space 


represents 


Che  spontaneous  P.O.P.  emission.  It  is  inde¬ 
pendent  of  the  temperature.  The  broadening 
factor  of  the  D.F.  l«  proportlen.il  to  distance 

and  Is  given  by  ~  t  V.”  It  is  an  in¬ 
creasing  function  pf  Che  temperature.  It  is 
seen  chat  alter  an  infinite  distance  from  the 
origin,  the  iJ.F.  reaches  its  equilibrium  profile. 
From  Eq.  I  the  .qull  iiiriun  D.F.  is  given  by 


P.o. 


il 
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We  present  a  novel  method  to  solve  explicitly  the  Boltzmann  equation  for  highly  energet¬ 
ic  electrons  interacting  with  polar  optical  phonons  and  scattering  mainly  in  the  forward 
direction.  In  this  approach,  the  collision  integral  of  the  Boltzmann  equation  is  reduced  to  a 
differential  operator  which  is  much  easier  to  manipulate  than  the  integral  form  and  does 
not  require  a  relaxation- time  approximation.  The  relaxation  of  the  distribution  function 
with  time  as  well  as  the  spatial  evolution  of  highly  energetic  electrons  are  calculated  and 
closed- form  expressions  for  the  distribution  function  are  given.  In  both  cases  the  behavior 
of  the  electron  distribution  is  characterized  by  two  fundamental  parameters:  a  drift  factor 
which  represents  the  net  rate  of  phonon  emission,  and  a  broadening  factor  which  is  propor¬ 
tional  to  the  latter  and  also  to  time  and  distance. 


I.  INTRODUCTION 

The  theory  of  high-field  transport  in  semicon¬ 
ductors  is  closely  related  to  the  solution  of  the 
Boltzmann  equation  for  high  carrier  energies.1  Ow¬ 
ing  to  its  complexity  (integro-differential  equation), 
the  Boltzmann  equation  cannot  be  solved  explicitly. 
In  the  post  two  general  methods  of  approximation 
have  been  proposed.  The  first  relies  on  the  concept 
of  electron  temperature  and  assumes  a  Maxwellian 
form  of  the  isotropic  part  of  the  electron  distribu¬ 
tion  function.2,3  Unfortunately,  this  analytical 
method  applies  only  to  the  case  of  very  high  elec¬ 
tron  densities.  The  other  methods  use  numerical 
techniques  (iterative  and  Moote  Carlo)  to  solve  the 
Boltzmann  equation.4,3  These  methods  are  more 
exact  but  are  rather  time  consuming  and  costly  and 
therefore  not  easily  applicable  to  semiconductor  de¬ 
vice  models. 

Recently,  transient  transport  phenomena  at  high 
energies  have  been  the  subject  of  considerable  in¬ 
terest  in  connection  with  “ballistic  transport’’  in 
very  short  devices.6  The  idea  is  that  for  small  de¬ 
vice  dimensions  (of  the  order  of  the  mean  free  path) 
the  charge  carriers  suffer  only  a  few  collisions  and 
gain  extremely  high  speeds.  In  the  calculations 
semiempirical  (Newton)  equations  such  as  the  fol¬ 
lowing  are  often  used7: 


d  e  mv* 
—mvd=qE — ~  , 
dt  rm 


(la) 


d  „  f-fo 

—e=qvdE - ,  (lb) 

dt  r, 

where  e  and  e0  are  the  average  and  the  zero-field 
energy,  respectively,  vd  is  the  drift  velocity,  E  the 
electric  field,  m  the  effective  mass,  and  rm  and  r, 
are  the  empirical  momentum  and  energy  relaxation 
times. 

This  method  has  been  subject  to  controversy  since 
the  boundary  conditions  have  been  oversimplified 
and  spatial  inhomogeneities  and  the  statistical  na¬ 
ture  of  the  charge  transport  have  been  neglect¬ 
ed.*'  |°  Moreover,  the  criteria  that  define  the  mean 
free  path  are  often  based  on  the  low-field  and 
steady-state  values  of  the  physical  parameters, 
whereas  the  calculations  are  applied  to  high-field 
and  transient  phenomena.  The  numerical  methods, 
Monte  Carlo,  etc.,  also  have  their  limitations.  In 
addition  to  their  high  cost,  they  cannot  easily  be  ap¬ 
plied  to  complicated  device  structures. 

In  this  paper  we  present  a  new  derivation  of  the 
Boltzmann  equation  for  fast  electrons  scattered  by 
polar  optical  phonons  (POP).  This  approach  can  be 
used  to  obtain  closed-form  integrations  of  the 
Boltzmann  equation  even  for  sophisticated  device 
structures,  provided  the  electrons  are  injected  at 
high  energies  and  the  electric  fields  away  from  the 
injecting  barrier  can  be  treated  as  perturbation.  In¬ 
jection  of  electrons  over  barriers  was  proposed  re¬ 
cently"  to  achieve  extremely  high  electron  veloci¬ 
ties  over  large  distances.  It  is  important  in  many 
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A  novel  method  previously  presented  to  explicitly  solve  the  Boltzmann 
equation  for  highly  energetic  electrons  interacting  with  polar  optical 
phonons  u  improved  in  order  to  take  into  account  the  non-parabolicity  of 
the  conduction  band.  For  the  case  of  field  free  spatial  diffusion,  the  carrier 
velocity  is  compared  with  Monte  Carlo  simulations.  For  GaAs  the  agree¬ 
ment  between  the  analytical  and  the  numerical  method  is  quite  good  at 
77  and  300  K.  Also,  the  results  derived  for  a  parabolic  band  structure  show 
a  strong  overestimation  of  the  carrier  velocity. 


IN  A  PREVIOUS  PAPER  [1  J,  we  presented  a  novel 
method  to  explicitly  solve  the  Boltzmann  equation  for 
highly  energetic  electrons  interacting  with  polar  optical 
phonons  (POP).  This  approach  is  based  on  the  consider¬ 
ation  of  forward  scattering  -  i.e.,  with  conservation  of 
the  electron  momentum  direction  -  due  to  the  electron- 
POP  interaction  at  high  energy.  The  collision  integral  of 
the  Boltzmann  equation  is  reduced  to  a  differential 
operator  which  is  much  easier  to  manipulate  than  the 
integral  form,  and  does  not  require  a  relaxation  time 
approximation. 

In  the  past  similar  methods  have  been  applied  for 
both  polar  [2, 3]  and  non-polar  optical  scattering  [4, 

5].  These  phenomenological  approaches  are  restricted  to 
energy  space  and  the  distribution  function  obtained  by 
these  methods  is  only  a  function  of  the  carrier  energy,  * 
not  of  the  carrier  momentum.  Our  present  method  is 
more  general  because  the  differential  operator  for  the 
POP  scattering  is  directly  derived  from  the  Boltzmann 
equation  at  high  energy.  Hence  we  do  not  introduce  any 
assumption  on  the  form  of  the  distribution  function 
which  contains  more  information  on  the  transport 
properties  (the  streaming  effect  for  example)  than  a 
distribution  defined  in  energy  space  only. 

The  differential  nature  of  the  collision  operator, 
which  is  similar  to  a  Fokker-Planck  equation,  is  well 
adapted  to  treat  the  stochastic  processes  of  non¬ 
equilibrium  phenomena  (6).  In  this  respect,  this  method 
appears  to  be  a  better  approximation  than  previous 
theories  which  exclude  statistical  considerations  as  well 
as  realistic  boundary  conditions  (ballistic  transport 
theory)  (7).  On  the  other  hand,  it  is  less  costly  than 
numerical  simulations  and  therefore  is  easily  applicable 
to  semiconductor  derice  models. 

The  basic  derivation  has  been  presented  before  [  I  ] 
assuming  a  parabolic  band.  The  importance  of  including 
the  nonparalolicity  for  high  energy  processes  is  obvious. 


We  restrict  our  considerations  to  scattering  pro¬ 
cesses  in  the  central  T-valley  and  assume  that  the  con¬ 
duction  band  is  isotropic,  omitting  any  consideration 
of  the  satellite  valleys.  This  approximation  is  reasonable 
for  energies  below  approximately  0.3  eV  in  GaAs. 

The  non-parabolicity  of  the  band  is  taken  into 
account  through  the  relationship  between  the  electron 
wavevector  k  and  its  energy  e 

— —  *  eil  +x«)  (1] 

1m 

where  in  agreement  with  Ehrenreich's  theory  [8) 


ef\  m0J 


Here  et  is  the  energy  gap  and  m  and  m0  are  the  effective 
mass  and  the  free  electron  mass,  respectively.  In  non¬ 
degenerate  statistics,  the  collision  integral  for  electron- 
POP  interaction  takes  the  form 
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where 


Here  o  is  Frdhlich’s  electron  -phonon  coupling  con¬ 
stant  and  is  the  POP  frequency  [*)J.  in  the  forward 
scattering  approximation  (I  J,  we  obtain 
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Investigation  of  Transient  Electronic  Transport 
in  GaAs  Following  High  Energy  Injection 

JEFFREY  YUH-FONG  TANG  and  KARL  HESS,  member,  ieee 


Atoowcf-W*  pnaeat  Mom*  Carlo  ranaiitio—  of  the  tnriwt  be¬ 
havior  of  dectroa*  injected  into  GaAs  at  Ugh  rairgira  and  accrtoratad 
(decelerated)  By  coaataat  dec  trie  fields.  Oar  cMcatefima  differ  f  toot 
pradooe  calcalatioaa  da*  to  the  indodon  of  a  raeiiitir  bead  stnetar* 
(oopiricd  paeadORotaatW)  atr  Ike  iejactioii  of  dectioo*  at  M^i 
•Mtjpra  (*4-.  via  a  heteiobairier).  The  raaaits  (hoar  that  a  earroar 
"cdBdon  ft**  window "  (CFW)  exists  with  respect  to  penwietsw  tech 
as  the  extsraal  electric  field,  the  ipjectioa  eaergy,  the  external  vdtagra, 
aad  the  seasicoadactor  dhaeashw.  WhMa  this  wiadow  avenge  dec- 
Boa  velocities  of  ~8  x  107  cra/s  caa  be  achieved  over  distaacas  of 
MT*  cm  ia  eautter  (soarce)-  aad  bear  Hire  stractaraa.  Voltaps  (field) 
parsiaeters  typical  for  the  collector  (dram)  ere  far  outside  the  CFW  sad 
allow  only  for  orach  isdaced  (by  a  factor  of  ~10)  elec  trow  velocities. 
We  also  discuss  thermal  noise  in  "balBstic  devices’*  aad  show  that  the 
noise  e<i Bivalent  temperatase  caa  be  exceedingly  high- 


I.  Introduction 

SINCE  THE  FIRST  calculation  of  transient  electron  trans¬ 
port  in  short-channel  FET’s  by  Ruch  [1],  overshoot  phe¬ 
nomena  of  the  electron  velocity  on  very  short  time  and  length 
scales  have  attracted  considerable  attention  [2 j  — {5] .  Shur 
[3]  and  Shur  and  Eastman  (5]  added  features  of  space  charge 
limited  conditions  to  the  velocity  overshoot.  They  investigated 
the  initial  transient  and  called  it  the  “near  ballistic”  regime. 
The  term  “ballistic”  is  difficult  to  define  and  is  currently 
applied  to  a  wide  range  of  device  parameters  and  dimensions. 
Investigations  in  (3]  and  [S]  indicate  that  “near  ballistic” 
transport  over  larger  distances  (> 1 000  A),  is  achievable,  necessi¬ 
tates  injection  of  electrons  at  higher  energy.  High  energy 
injection  was  also  discussed  by  Hess  [6]  for  low  temperature 
transport  free  of  scattering  events  over  extremely  large  distances 
L  >  10~4  cm  (for  electron  energies  below  0.036  eV). 

The  criteria  for  designing  devices  and  choosing  materials  such 
that  high  transient  speeds  can  be  advantageously  achieved  using 
high  energy  injection  are  the  key  issue  of  this  paper.  We 
approach  the  problem  by  a  Monte  Carlo  simulation  which 
includes  the  details  of  the  band  structure  as  calculated  by  the 
empirical  pseudopotential  method.  The  electron  is  started  at 
higher  energy  values  (not  at  the  bottom  of  the  band  as  done 
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by  Ruch)  to  simulate  injection  over  a  (hetero)barrier.  The 
electric  field  is  chosen  to  be  constant  over  the  whole  distance 
of  the  simulation,  in  contrast  to  the  choice  of  the  self-consis¬ 
tent  field  by  Shur  (3] .  The  reason  for  this  is  the  following: 
Within  the  CFW,  the  electron  velocity  does  not  vary  strongly 
over  the  device  length  because  the  electrons  are  already  injected 
at  high  velocities  (in  contrast  to  the  cases  considered  by  Shur). 
Therefore,  the  carrier  concentration  is  rather  constant  and  the 
electric  field  induced  by  the  carriers  is  of  minor  importance. 
In  all  other  respects  our  calculation  is  standard.  We  do  not 
include  the  intracollisional  field  effect  [7]  because  it  would 
greatly  complicate  the  computations.  In  the  CFW  it  also 
would  make  only  minor  contributions  since  the  electric  fields 
in  this  window  are  small. 

With  respect  to  device  applications,  the  result  of  the  calcula¬ 
tions  can  be  summarized  as  follows:  On  a  length  scale  of  1000 
A,  emitter  (source)-  and  base-like  structures  may  show  effects 
typical  for  collision-free  transport  (if  carefully  designed); 
coilector(drain>iike  structures  will  not. 

n.  Physical  Model  and  Method  of  Computation 
As  discussed  above,  we  consider  high  energy  electrons  injected 
into  GaAs  (e.g.,  from  Al^Ga^As  or  6 -like  electric  fields 
created  by  spice-charge  layers).  The  transition  is  assumed  to 
be  abrupt,  i.e.,  the  electron  gains  the  kinetic  energy  A E  and  a 
forward  momentum  AS  when  transferring  to  the  GaAs  with¬ 
out  any  energy  losses. 

As  illustrated  below,  the  calculations  of  the  self-consistent 
field  resulting  from  external  voltages  and  electron  redistribu¬ 
tion  in  tha  GaAs  is  a  difficult  problem  and  depends  on  many 
details,  most  importantly: 

1)  the  boundary  conditions  of  the  injection, 
ii)  the  statistics  of  electrons  and  impurities  in  the  GaAs, 

Ui)  the  level  of  injection,  and 

hr)  the  velocity  distribution  and  the  injection  energy. 

Let  us  assume  that  our  device  is  short  in  the  * -direction  but 
rather  wide  in  they -direct  ion.  We  then  can  still  define  average 
quantities  such  as  the  density  of  electrons  in  a  meaningful  way 
and  use  a  continuum  picture  as  follows:  If  the  distribution 
function  is  denoted  by /we  define  electron  density 


2 


f»f 


or  the  current  density  as 
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The  dynamics  of  transient  transport  and  transferred  electron  behavior 
are  studied  under  the  condition  of  high-energy  electron  injection.  This 
study  makes  use  of  a  Monte  Carlo  simulation  with  the  unique  inclusion  of 
realistic  band-strurture  as  derived  from  an  empirical  pseudopotential  method. 
The  r-L-X  valleys  are  implicitly  coupled,  thereby  permitting  the  study  of 
transferred  electron  behavior  in  a  natural  way.  The  details  of  the  r-L-X 
intervalley  electron  transfer  characteristics  are  presented  showing,  for  the 
first  time,  the  dynamical  interplay  between  the  X  and  L  valleys,  it  is 
readily  seen  from  our  results  that  electron  transport  in  gallium  arsenide 
reauires  a  multivalley  description;  in  addition,  the  onset  of  inter-valley 
electron  transfer  Is  snown  to  depend  strongly  upon  the  electron  injection 
energy. 


Gallium  arsenide  electronic  devices  with 
iransit  lengths  in  the  submicron  and  ultrasub- 
“icron  regions  are  coming  to  fruition  due  to  the 
advent  of  MBEW  and  M0CVD.3  The  Intention  of 
:nis  study  Is  to  describe  a  variety  of  transient 
transoort  characteristics  pertinent  to  gallium 
arsenide,  such  as  high-speed  transport  resulting 
from  high-energy  electron  Injection,  and  the  dy¬ 
namics  of  transferred  electron  behavior. 

The  transport  properties  described  herein 
are  obtained  by  means  of  a  Monte  Carlo  simula¬ 
tion.  The  Monte  Carlo  method  utilized  In  this 
;aoer  is  unique  In  that  It  Includes  a  realistic 
*and-structure  as  derived  from  an  empirical 
pseudopotential  method.  In  this  simulation,  the 
*-l-X  valleys  are  implicitly  coupled,  thereby 
Permitting  the  study  of  transferred  electron  ef¬ 
fects  in  a  natural  way. 

It  is  found  that  the  transient  velocities 
achievable  In  gallium  arsenide  can  be  quite 
•r,*e  (3-10  X  107  cm/sec);  however,  achievement 
::  Mjc.n  high  velocities  is  shown  to  depend  sen- 
-:*.ivniy  upon  the  electric  field,  the  injection 
■''•my,  and  the  transit  distance.  For  devices, 
they  are  currently  understood,  we  show  that 
"  'reeds  can  be  achieved  only  in  base- 
'iirce  emitter-)  like  structures;  it  is  noted 
" ,r  in  collector-  (drain-)  like  structures  the 
; ■_••• r.-rwsent  hioh-voltage  drops  do  not  permit 
•■■ml Testation  of  high-velocity  (above  satu- 
: ion. velocity)  transport. 

’Iso,  >/e  discuss  In  detail  the  r-X-l  inter- 
i11«v  -'loctron  transfer  characteristics  showina, 
1  '*e  first  tire,  the  dynamical  interplay  be¬ 


tween  the  X  and  L  valleys.  It  is  readily  seen 
from  our  results  that  electron  transport  in  gal¬ 
lium  arsenide  generally  requires  a  multivalley 
description;  In  addition,  the  intervalley  elec¬ 
tron  transfer  Is  also  shown  to  be  strongly  depen¬ 
dent  upon  the  electron  Injection  energy. 

The  transport  properties  described  herein 
were  obtained  by  means  of  a  Monte  Carlo  simula¬ 
tion.  The  details  of  the  calculations,  to  with¬ 
in  a  few  minor  adjustments,  have  been  discussed 
elsewhere;4 *5 therefore,  an  In-depth  discussion 
of  the  simulation  is  omitted  here.  In  general, 
the  calculation  makes  use  of  the  following  as¬ 
sumptions: 

(1)  The  electrons  propagate  between  scat¬ 
tering  events  in  accordance  with  the  equations 
of  motion 

hfc  *  eE  (1) 


•  (2) 

«*• 

Here  E,  the  applied  electric  field,  is  constant 
in  space  and  time,  and  is  assumed  to  be  turned 
on  instantaneously;  also,  c(k)  is  the  energy- 
band  dispersion  relation  calculated  by  the  em¬ 
pirical  pseudopotential  method,'  but  with  proo- 
erly  adjusted  r-L-X  intervalley  spaclngs  (0.33eV 
and  O.S 2eV,  respectively). 

(11)  Folar  optical  scattering,  acoustic  de¬ 
formation  potential  scattering,  and  intervalley 
scattering  mechanisms  are  taken  into  account; 
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ABSTRACT 

Traaaioat  traaaport  of  alaetroaa  ia  ZaP  is  stadiad  nadar  tba  condition  of 
higir-eaergy  alaetroa  iajaetioa.  Tbit  study  aakas  usa  of  a  Monte  Carlo 
aiamlatioa  with  tba  caique  iaolaaioa  of  realistic  baad~structnre  as  derived 
tram  aa  empirical  pseudopotential  method.  The  results  obtained  herein  for  InP 
are  qualitatively  siailar  to  those  previously  obtained  by  the  authors  for 
OaAs.  Quantitatively,  it  is  found  that  ultrahigh  alaetroa  drift  velocities 
( a  10*  ca/sac)  persist  for  nueh  higher  electric  fields  and  over  aueh  longer 
distance  of  electron  traversal  ia  Za?  as  compared  to  GaAs. 
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ABSTRACT 

fe  present  a  Monte  Carlo  simulation  of  the  behavior  of  a  collecting  n-n+ 
contact  in  IaP.  Electrons  are  injected  at  high  energy  into  the  In?  and  are 
accelerated  by  an  applied  electric  field  over  a  length  of  *  1000  X.  At  the 
collecting  contact  they  encounter  a  possible  reflection  bach  into  the  device. 
The  reflection  coefficient  at  the  contact  is  chosen  to  vary  between  0  and 
0.70.  The  results  show  that  the  average  electron  drift  velocity  is  greatly 
lowered  at  low  fields  throughout  the  entire  device  by  the  reflection  at  the 
contact. 
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ABSTRACT 

Caloalations  of  tba  steady  state  and  transient  electron  drift  velocities 
and  iapaet  ionization  rate  are  presented  for  GaAs*  InP,  and  InAs  based  on  a 
Monte  Carlo  simulation  using  a  realistic  band  structure  derived  from  an 
eapirical  pseudopotential.  The  inpact  ionisation  results  are  obtained  using 
collision  broadening  of  the  initial  state  and  are  found  to  fit  the 
experimental  data  well  through  a  vide  range  of  applied  fields.  In  InP  the 
iapaet  ionisation  rate  is  aueh  lover  than  in  GaAs  and  no  appreciable 
anisotropy  has  been  observed.  This  is  due  in  part  to  the  larger  density  of 
states  in  InP  and  the  corresponding  higher  electron-phonon  scattering  rate. 

The  transient  drift  velocities  are  calculated  under  the  condition  of  high 
energy  injection.  The  results  for  InP  shov  that  higher  velocities  can  be 
obtained  over  1000  -  1500  1  device  lengths  for  a  aueh  larger  range  of 
launching  energies  and  applied  electric  fields  than  in  GaAs.  For  the  ease  of 
InAs*  due  to  the  large  iapaet  ionisation  rate*  high  drift  velocities  can  be 
obtained  since  the  ionisation  acts  to  liait  the  transfer  of  electrons  to  the 
satellite  ainiaa.  In  'Che  absence  of  iapaet  ionisation,  the  electrons  shov  the 
usual  runavsy  effect  and  transfer  readily  occurs*  thus  lovering  the  drift 
velocity  substantially. 


Theory  of  hot  electron  er lesion  from  silicon  into  silicon  dioxide 

J.  Y.  Tang*1  and  Karl  Hess 

Department  pf Electrical  Engineering  and  Coordinated  Science  Laboratory,  University  of  Illinois  at  Urbana- 
Champaign,  Urbana,  Illinois  61301 

(Received  4  April  1983;  accepted  for  publication  1  June  1983) 

We  present  Monte  Carlo  simulations  of  the  hot  electron  emission  from  silicon  into  the  oxide  of 
metal  oxide  silicon  transistors.  The  calculations  include  the  pseudopotential  band  structure  and 
quantum  effects  such  as  collision  broadening  due  to  the  electron— phonon  interaction.  As  a  result, 
we  present  a  set  of  transport  parameters  which  well  describes  all  hot  electron  effects  in 
(including  saturation  velocity  and  impact  ionization).  We  also  show  that  the  «sliiairt«.  broadening 
effect  leads  to  an  effective  barrier  lowering  and  may  require  that  voltages  be  trahwl  down  fig  below 
the  interface  barrier  height,  of  —>3.1  V  in  order  to  avoid  hot  electron  '****i— *'**' 

PACS  numbers:  72.20.Ht,  72.10.Bg,  85.30.Tv,  73.40.Qv 


I.  INTRODUCTION 

Electronic  transport  in  the  Si/Si02  system  has  been 
studied  extensively  because  of  its  important  role  in  silicon 
device  technology.  One  of  the  most  intriguing  problems  of 
hot  carrier  transport  at  the  Si/Si02  interface  has  been  the 
emission  of  hot  electrons  or  holes  from  silicon  into  SiO,.1"4 
Various  experiments  have  been  performed  to  study  this  ef¬ 
fect.  From  our  viewpoint,  Ning’s  experimental  setup  is  the 
most  interesting.  Ning  et  al?A  measured  the  absolute  emis¬ 
sion  probability  of  electrons  which  were  optically  generated 
in  the  silicon  depletion  layer  and  accelerated  toward  the  Si- 
SiO:  interface.  Past  models  explaining  the  emission  process 
involved  many  simplifying  assumptions,  for  example,  an  en¬ 
ergy  independent  mean  free  path  and  a  parabolic  band  struc¬ 
ture.  At  energies  of  3.0  eV  above  the  conduction  band  edge, 
the  band  structure  is,  of  course,  highly  nonpafabolic.  Trans¬ 
port  at  energies  at  which  the  emission  takes  place  must  also 
consider  different  valley  types  {XJ. )  and  more  than  one  con¬ 
duction  band.7  Moreover,  impact  ionization  definitely  plays 
an  important  role. 

It  is  the  purpose  of  this  paper  to  present  a  Monte  Carlo 
simulation  of  the  hot  election  emission  process  which  in¬ 
cludes  a  realistic  band  structure  and  quantum  effects  such  as 
collision  broadening.  It  will  be  shown  that  the  collision  broa¬ 
dening  effect  introduces  an  effective  barrier  lowering  which 
is  important  for  considerations  of  “scaling  away”  the  hot 
electron  effect 

IL  SUMMARY  OP  EXPERIMENTAL  RESULTS 

The  devices  used  by  Ning  et  al?  were  n -channel  polysi- 
lieon-SiOj-Si  field-effect  transistors.  Electrons  are  optically 
generated  and  injected  into  Si02  as  illustrated  in  Figs.  1(a) 
and  1(b).*  Source  and  drain  were  grounded,  a  negative  bias 
was  applied  to  the  substrate,  and  a  positive  bias  was  applied 
to  the  gate.  The  advantage  of  this  arrangement  is  that  the 
gate  voltage  and  the  substrate  voltage  can  be  varied  indepen¬ 
dently.  OpticaUygenerated  electrons  which  diffuse  into  the 
depletion  region  are  accelerated  toward  the  Si-Si02  inter¬ 


face.  The  majority  of  camera  that  do  not  overcome  the  inter¬ 
face  barrier  are  collected  aad  give  riec  to  source  and  drain 
currents  ID.  The  camera  that  overcome  die  barrier  contri¬ 
bute  to  the  gate  current  Ia  if  they  are  not  trapped  in  the 
silicon  dioxide  layer.  By  measuring  the  gate,  the  source,  and 
the  drain  currents,  the  absolute  emission  probability  can  be 
obtained  from 


where  7toul  —2 ID  is  the  total  current  from  the  substrate. 

The  doping  profile  in  the  silicon  substrate  can  be  ap¬ 
proximated  by  a  Gaussian  distribution 

Na  (■*)  =  JV*  +  C0  exp^  —  ,  (2) 

where  x  is  the  distance  from  the  Si-Si02  interface,  N„  is  the 
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HO.  t.  Experimental  arrangement  (a)  and  schematic  band-structure  dia¬ 
gram  |b|  defining  the  hot  electron  emission  from  silicon  into  silicon  dioxide. 
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ABSTRACT 

Ve  preaant  a  Mont a  Carlo  simulation  of  alaetronie  transport  in  new  forma 
of  GaAs-OaAlAs  batarostrnoturas  dasignad  for  achieving  high  drift  velocities. 
Average  alaetron  speeds  of  ~  5.0x10^  ea/sao  tbrongb  an  antira  structure  length 
of  0.5  |ia  have  been  calonlsted.  This  represents  a  marked  improvement  in  spaed 
over  a  single  barrier  emitter  atrsetnre.  The  basie  physical  principle  is  the 
following:  the  electrons  are  confined  to  the  gamma  valley  by  losing  excess 
kinetic  energy  gained  from  an  overlaid  accelerating  field.  The  mechanism  for 
the  energy  loss  is  a  series  of  ascending  potential  steps.  In  this  way 
transfer  to  the  subsidiary  minima  is  avoided  and  very  high  velocities  are 
possible.  Paradoxically,  this  means  that  because  the  electrons  lose  kinetic 
energy  their  velocity  remains  high  or  actually  increases. 


scheme.  At  was  noted  in  [6|,  theti  EP  calculation!  (4| 
differ  from  the  experimental  lerdtt  by  about  60meV 
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